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Abstract 



We derive exact Friedmann-Robertson- Walker cosmological solutions in gen- 
eral scalar-tensor gravity theories, including Brans-Dicke gravity, for stiff 
matter or radiation. These correspond to the long or short wavelength modes 
respectively of massless scalar fields. If present, the long wavelength modes of 
such fields would be expected to dominate the energy density of the universe 
at early times and thus these models provide an insight into the classical 
behaviour of these scalar-tensor cosmologies near an initial singularity, or 
bounce. The particularly simple exact solutions also provide a useful exam- 
ple of the possible evolution of the Brans-Dicke (or dilaton) field, eft, and the 
Brans-Dicke parameter, to{4>), at late times in spatially curved as well as flat 
universes. We also discuss the corresponding solutions in the conformally 
related Einstein metric. 
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I. INTRODUCTION 



If we hope to describe gravitational interactions at energy densities approaching the 
Planck scale it seems likely that we will need to consider Lagrangians extended beyond the 
Einstein-Hilbert action of general relativity. The low-energy effective action in string theory, 
for instance, involves a dilaton field coupled to the Ricci curvature tensor Scalar fields 
coupled directly to the curvature appear in all dimensionally reduced gravity theories, and 
their influence on cosmological models was first seriously considered by Jordan M. These 
models have been termed scalar-tensor gravity, the best known of these being the Brans- 
Dicke theory ||. Gravity Lagrangians including terms of higher order in the Ricci scalar 
can also be cast as scalar-tensor theories with appropriate scalar potentials. 

The belief that modified gravity theories may have played a crucial role during the early 
universe has recently been rekindled by extended inflation ||. In this scenario a scalar- 
tensor gravity theory allows the first order phase transition of the "old" inflationary model 
|7]] to complete. This arises because the scalar field, 0, (henceforth the Brans-Dicke field, 
essentially the inverse of the Newton's gravitational "constant") damps the rate of expansion 
and, in the original extended inflationary scenario based on the Brans-Dicke theory, turns 
the exponential expansion found in general relativity into power law inflation ||. However, 
Brans-Dicke theory is unable to meet the simultaneous and disparate requirements placed 
by the post-Newtonian solar system tests || and by the need to keep the sizes of the bubbles 
nucleated during inflation within the limits permitted by the anisotropics of the microwave 
background flOH . 

This situation may be averted through the consideration of more general scalar-tensor 
theories in which the parameter u, a constant in Brans-Dicke theory, is allowed to vary as 



uj((j)) JTI|. But in such cases one needs to understand better the cosmological behaviour of 
these general theories, in order to assess their implications on our models of the universe. 
Direct observations mainly constrain these theories at the present day in our solar system 



T2|| imposing a lower bound u > 500 and requiring that u 3 (du / d<p) should approach zero. 



On a cosmological scale, the principal limits arise from the consideration of effects upon the 
synthesis of light elements, indicating that at the time of nucleosynthesis similar bounds 
hold Jl3| . 



Most of the work that can be found in the literature on solutions of scalar-tensor theories 
concerns the particular case of Brans-Dicke theory fl4Hl9|. The properties of more general 
scalar-tensor cosmologies have been discussed recently and exact solutions derived 

for the vacuum and radiation models JZIJ (where p = p/3) corresponding to the particular 
situation where the scalar field is sourceless (because the matter energy-momentum tensor 
is traceless). In this paper we show how to extend this method to derive exact solutions for 
the homogeneous and isotropic cosmological models with a perfect fluid characterized by the 
equation of state p = p, which does act as a source for the Brans-Dicke field. These models 



represent the evolution a homogeneous massless scalar field ||23|| . Such a scalar field may 
describe the evolution of effectively massless fields, including in the context of superstring 
cosmology the antisymmetric tensor field which appears in the low energy string effective 
action gH. 

As the energy density of a barotropic perfect fluid with p = (7 — l)p evolves as p oc a -37 , 
such "stiff matter" would be expected to dominate at early times in the universe [p5f over 



2 



short wavelength modes or any other matter with p < p. Thus our solutions provide an 
important indication of the possible early evolution of scalar-tensor cosmologies. 

As in , our solutions will be given in closed form in terms of an integration depending 
on c<j(0), which can be performed exactly in many cases and numerically in all cases. The 
general field equations are given in Section for scalar-tensor gravity and we solve these in 
Friedmann-Robertson- Walker metrics for vacuum, stiff fluid and radiation models in Sec- 
tion JTJ. The equivalent picture in the conformally transformed Einstein frame is presented 
in Section [IV|. Conclusions on the general behaviour of solutions are presented in our final 
section. 



II. SCALAR TENSOR GRAVITY THEORIES 



The scalar-tensor field equations |26| are derived from the actionQ 

S = / C^Xyf^g 

167T Jm 



4>{R- 2A(0)) - ^p-g^afo + 16<K£ m 



(2.1) 



where R is the usual Ricci curvature scalar of the spacetime, <fi is the Brans-Dicke scalar 
field, u>((f)) is a dimensionless coupling function and C m represents the Lagrangian for the 
matter fields. It is clear that the scalar field plays the role which in general relativity is 
played by the gravitational constant, but with now a dynamical variable. 

The particular case of Brans-Dicke gravity arises when we take uo to be a constant and 
A = in the Lagrangian of Eq. ( j2.1|) . The A(0) potential is the natural generalisation of 
the cosmological constant A. It introduces terms which violate Newtonian gravity at some 
length scale. In what follows we will leave u)(<f>) as a free function but consider only models 
in which A is zero. This should be valid at least at sufficiently early times when we expect 
kinetic terms to dominate, and will also avoid introducing too many free functions into our 
analysis. (Note that the Lagrangian is sometimes written in terms of a scalar field ip with a 
canonical kinetic term so that <fi = f(<p) and cj(<f>) = f/2(df /dip) 2 .) 

Taking the variational derivatives of the action ( |2.1| ) with respect to the two dynamical 
variables g a b and (ft and setting A (</>) = yields the field equations 



Rab — -^9abR 



□c 



~6 V a b ~ 2 9ab9 ! ' 



+- (V a V b - g ab u 
I r„ ~ c rf 



2a;(0) + 3 

where T = T° is the trace of the energy-momentum tensor of the matter defined as 



(2.2) 
(2.3) 



One can also include an additional boundary term dependent on the extrinsic curvature of the 
boundary M, as is required in general relativity [^7| to allow for the variation of g a b,c ° n the 
boundary. 
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ab 



(2.4) 



It is important to notice that the usual relation VfeT ab = establishing the conservation 
laws satisfied by the matter fields holds true. This follows from the assumption that all 
matter fields are minimally coupled to the metric g a t which means that the principle of 
equivalence is guaranteed. The role of the scalar field is then that of determining the 
spacetime curvature (associated with the metric) produced by the matter. Matter may be 
a source of the Brans-Dicke field, but the latter acts back on the matter only through the 
metric [EM . 



III. FRIEDMANN ROBERTSON WALKER MODELS 

We consider homogeneous and isotropic universes with the metric given by the usual 
Friedmann-Robertson- Walker (FRW) line element 



ds 2 = -dt 2 + a{tf 



dr 2 



+ r 2 (d6 2 + sinO 2 dip 2 " 



1 — kr' 



(3.1) 



The field equations for a scalar-tensor theory, where we allow the coupling parameter u 
to depend on the scalar field 6, but restrict the potential A to be zero, are then 



6 6 S 1 a 2 3 6 



(3.2) 



+ 



3 <i u{4>) 
a 2uj(6) + 3 



8vrp 



2a; (6) + 3 



(4-37) 



(3.3) 



H + H 2 + ^^-H. 



l-Kp (37-2)^ + 3 1 



to 



2uj{6) + 3 2 2u{6) + 3 



(3.4) 



These equations differ from the corresponding equations of Brans-Dicke theory through the 
presence of terms involving Co in the two latter equations. 

Several authors have studied cosmological solutions of the Brans-Dicke theory for a FRW 
universes filled with a perfect fluid P,|T4}-|T7|] . Nariai |2S| derived power law solutions for the 
flat FRW universe with a perfect fluid satisfying the barotropic equation of state p = (7 — l)p, 
with 7 a constant taking values in the interval < 7 < 2. 

The solutions of these equations of motion are defined by four integration constants 
whereas the corresponding solutions in general relativity depend on only three ||18|| . In fact, 
in addition to the values of a(t ), a (to) an d 4>(t ) oc 1/G(t ) , we now need p(t ) (or <fi(t ) 
instead) as well. Originally, as done by Brans and Dicke, and by Nariai, this extra freedom 
was eliminated by requiring that 0a 3 should vanish when a approaches the initial singularity 
at a = 0. In a flat FRW model this restricts one to obtaining only the power law solutions 
of the Brans-Dicke theory. Here we shall keep our analysis more general. 
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The derivation of the general barotropic Brans-Dicke solutions for the spatially flat 
(k = 0) model was done by Gurevich, Finkelstein and Ruban [16]. Their solutions, which 
cover all the space of parameters of the theory, render clear a very important feature of the 
behaviour of the flat cosmological models, namely that for the solutions which exhibit an 
initial singularity (those with u > —3/2 ) the scalar field dominates the expansion at early 
times, whilst the later stages are matter dominated and approach the behaviour of Nariai's 
solutions for uj > 2(7 — 5/3)/(2 — r y) 2 . Note therefore that the stiff fluid solutions are unique 
in that they do not approach Nariai's power law solutions at late times other than in the 
limit uj —>■ 00. 

The same solutions for the flat model in the cases of vacuum, stiff matter (p = p) and 
radiation (p = p/3) were rederived later by Lorentz-Petzold [17j using a different method 
which enabled him to also obtain solutions for the non-flat models. We use this method in 



an improved form |22| to derive solutions for the general scalar-tensor theories. 
Using the conformal time variable r\ defined by the differential relation 

dt = a drj , 



(3.5) 



and the variables 



X = 



(3.6) 



and 



Y 



'2uj + 3 



(3.7) 



we can re-write the above field equations as 

(X') 2 + 4 k X 2 - (V X f = AM X a 4 " 37 



(3.8) 



[Y'X]' = M(4 - 37) a 4 " 37 , (3.9) 



X" + 4 k X = 3 M(2 - 7) a 4-37 , (3.10) 

where the density p = 3M/87ra 37 for a barotropic fluid with M a constant. The prime 
denotes differentiation with respect to 77. Our variables are akin to those used by Lorentz- 



Petzold [17] when solving for the Brans-Dicke theory. To that extent the method we explore 
here is a generalization of his method of obtaining decoupled equations. Note that whenever 
X is negative this must correspond to a negative value for 0. In what follows, unless 
otherwise explicitly stated, we shall assume that u > —3/2 to guarantee the positiveness 
of the function under the square root in Eq. ( |3.7| ), although it would be straightforward to 
redefine Y(<p) for the case of uj < — 3/2. 

This system considerably simplifies for the two particular cases: 7 = 4/3 (radiation) and 
7 = 2 (stiff matter). We shall show in the next section precisely how these correspond to 
the short and long wavelength limits of a massless field. In either case the full integration is 
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again possible, provided we specify the function uj{4>). For other values of 7, the equations 
carry an explicit dependence on a which cannot be integrated by the method adopted here. 
An alternative approach for these latter cases based on another method of integration of the 



original field equations is presented elsewhere [pO 



Anisotropic cosmologies have also been considered in the literature, again principally for 
Brans-Dicke gravity. We will show elsewhere how our method may be extended to derive 
solutions for general scalar-tensor gravity in anisotropic models piT . 



A. Scalar field evolution 

A minimally coupled scalar field, a, whose energy-momentum tensor 

T ab = (g c a gi - -g a bg cd ^j , 

= (p + p)u a u b + pg ab , (3.11) 
corresponds to a perfect fluid with density 

p = p = ^\g ab a a a b \, (3.12) 



and normalised velocity field 



\9 cd cr,c^,d\ ' 



The scalar field itself obeys the wave equation 

□a = , (3.14) 

which in a FRW metric reduces to 

3 

- a - 3H& + 9 ij V i ^3 a = ■ ( 3 - 15 ) 

If we consider plane wave solutions of the form a = <J q (t) exp(i J2 q%x l ) then for (q/a) 2 ^> 
H 2 , H, k/a 2 this can be re-written as 

(aa q )" + q 2 (aa q ) = 0, (3.16) 

where q 2 = J2q 2 - This corresponds to the usual flat space result for plane waves where 
(aa q ) cx exp(iqrj) and thus = q a q b cr 2 oc a -4 with the null four-vector q a = (q,qi). This 
is an anisotropic stiff fluid but if we consider an isotropic distribution of short wavelength 
modes averaged over all spatial directions this produces a perfect fluid with (p) = p/3, the 
usual result for isotropic radiation. 

For long wavelength modes in an FRW universe (with comoving wavenumber (q/a) 2 <^ 
H 2 , H) we can neglect spatial gradients in the field and the first integral of Eq. ( |3.16| ) yields 
a 3 a =constant and thus p = p = a 2 /2 oc a -6 , i.e. a stiff fluid. 
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Clearly the dividing line between these long and short wavelength modes changes as the 
comoving Hubble length or curvature scale evolves^]. In a conventional (non-inflationary) 
cosmology the comoving Hubble length shrinks as we consider earlier and earlier times in 
an expanding universe so that as a — > all modes must be "outside the horizon" and evolve 
as a homogeneous stiff fluid lending support to our contention that the stiff fluid solutions 
will be important in determining the classical behaviour of scalar-tensor cosmologies near 
any initial singularity. In any case, as already remarked, the energy density of a barotropic 
perfect fluid evolves as p oc a -37 and so the energy density of a stiff fluid will eventually 
dominate as a — > over any matter with a barotropic index 7 < 2. 

In what follows we shall consider only the extreme short and long wavelength modes of 
the massless field, neglecting the intermediate regimes. 



B. Vacuum solutions 



Let us first consider the field equations in vacuum; 

(X') 2 - (Y' X) 2 + 4k X 2 = , 
(Y' X)' = , 
X" + 4 k X = . 



(3.17) 
(3.18) 
(3.19) 



Both Eqs. ( |3.1S| ) & ( |3.19| ) are easily integrable, and X{j]) is independent of the particular 
uj{4>) dependence. 

Solving Eq. ( ff39p yields 



x(v) = 



A 7] 

4 sin (277) 
t ~ sinh (277) 



for k = , 
for A; = +1 
for k = — 1 



(3.20) 



with A an arbitrary integration constant [see Fig. In what follows we will find it most 
useful, and succinct, to write this as 



x( v ) 



±At 

1 + kT 2 



(3.21) 



in terms of the new time variable 

t(t]) = 



\V\ 
I tan 77 1 
tanhr/l 



for k = , 
for A; = +1 
forfc = — 1 



(3.22) 



2 Indeed this is precisely how long wavelength perturbations are produced in the inflaton field 
from originally short wavelength vacuum fiutuations as the comoving Hubble length shrinks during 
inflation J32| |. This highlights the potential importance of quantum effects which we shall neglect 
in this purely classical treatment. 
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It is convenient to define r as a non-negative quantity and choose the plus or minus sign 
in Eq. ( p. 21 ) according to whether rj is greater or less than zero respectively. This only 



amounts to a different choice of the integration constant and so can be absorbed in our 
choice of A. In practice, because a 2 = X/ <fi must always be non-negative only one choice of 
±A corresponds to a real solution anyway. For the allowed choice of A, r may then either 
increase or decrease with conformal (and thus also with proper) time. 
From Eq. (|3.18|) we obtain 



Y' X = f = const . (3.23) 



This latter result implies 



Now, notice that Eq. (|3.19|) has the first integral 

X' 2 + Ak X 2 = A , (3.25) 



where A = A 2 for the solutions given in Eq. (|3.21|) . Thus, substituting Eq. ( |3.23|) into 



Eq. fl3.17| ) we obtain a relation between the constants / and A such that A = ±/. 

Given the definition of X, and the fact that we know X(j]) from Eq. (|3.20|) , we realize 
that provided we know the particular form of uj{4>) we can obtain (j)(rj) from Eq. (|3.24f) , and 
then derive the scale factor a(rj) as 

a 2 (r]) = -X . (3.26) 



To obtain we have to invert Y ((/>), given by the left-hand side of Eq. ( |3.24 ) and use the 
fact that we know the right-hand side 

f f 

Y(rj) = — drj =± In t(t]) + constant. (3.27) 
J X 

We see that as r — > or r — > oo (at early or late times) the function Y must diverge. For 
instance in the case of Brans-Dicke gravity where u is a constant this implies that <fi — > 
or <p —>■ oo. However, there is a priori no prescription for uj(<f>). Thus, we are led to consider 
some specific u(4>) dependences which hopefully will shed some light onto general results 
concerning the dependence of the solutions on the form of a> (</>). 

Even without solving these equations for a particular u(<p) we can come to some general 
conclusions about how these vacuum solutions behave. As X — >• and (X'/X) 2 — ► oo the 
curvature becomes negligible in Eq. (|3.17| ) and we see that 




^3.28) 



Thus using the Eq. (|3.7| ) we have 
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(3.29) 
(3.30) 

and the initial singularity (with a — > ±00) can only be avoided for u — > 0. 

A necessary condition for any turning point in the evolution of the scale factor is 

Thus for k < a turning point can only occur when u < 0. This corresponds to the (V0) 2 
term in the action of Eq. ( |2.1| ) having the "wrong sign" , in that it can contribute a negative 



effective energy density. Turning points can occur in closed models even if u> > 0, just as 
they can occur in general relativity. Note that the sign of the gravitational coupling, 0, (and 
thus X) is irrelevant in this vacuum case. 



Vacuum solutions in Brans-Dicke gravity 

Let us consider first the case u(4>) = ujq =constant, corresponding to the Brans-Dicke 
theory. Then 



Y = J 2 -^±^\n± } (3.32) 



and thus 



±0 



■ 



(3.33) 

A r lT/3 

= ' (M4) 



where we have written (3 = y3/(2c<j + 3). These solutions are plotted in Fig. Q and Fig. |3|. 

The k = solutions correspond to those derived by O'Hanlon & Tupper If we 

convert them to proper time they read a(t) = aot q± and <fi = 9± \ where 

q± = —, — ; r , (3.35) 

3( W +1±V^) 



[see Fig. §]. The k 7^ solutions were obtained by Lorentz-Petzold Jl7] and by Barrow |p2 



As 77 — > 0, and thus a — > for > 0, they approach the k = power law behaviour. Also 
note that all solutions exhibit two branches. This is a consequence of the identity A = ±/ 
between the integration constants. Each branch corresponds to different signs of <p/4>. In 
fact, the q + branch is associated with an increasing |0|, which means that G approaches zero 
in the t — > 00 limit. Since this branch corresponds to a slower expansion, we shall follow 



Gurevich et al [16] in calling it the slow branch. On the contrary, the g_ fast branch has a 



decreasing \(f>\ and G, consequently, approaches ±00 with time. 
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Note that r — > (and thus 77 — > 0) coincides with a — > for both branches if and only if 
u; > (and thus (3 < 1), in agreement with our earlier arguments. For u; < the solutions 
do not have zero size at r = but are still singular in the sense that the Ricci curvature 
scalar, for instance, diverges. 

We can choose 0o to be either positive or negative and thus the sign of the gravitational 
"constant" is arbitrary as we would expect for solutions of the field equations in vacuum. Of 
course a 2 must remain positive so the product we require 4>qA > 0. Because of our definition 
of r in Eq. ( |3.22| ) we also have two distinct solutions corresponding to whether r decreases 
with time, corresponding to rj < and a collapsing universe as r — ► for u > 0, or increases 
with rj > for a universe expanding from r = if uj > 0. 



Vacuum solution with uj — » 00 
The simplest function which includes a divergent u;(0) at a finite value of = 0*, is 

2cu(0) + 3 = (2u + 3) . (3.36) 

<P* - <P 

The integral in Eq. (|3.7|) then yields 



Y(<f>) = l n ( v^" r J = ± lnr + constant (3 _ 37) 



which in turn gives 



(3.38) 



a 2 = — ± - R '- — , 3.39 



[see Fig. |5| where we have written (3 = ±y3/(2a;o + 3), although in fact the choice of ± is 
irrelevant here for 7^ 0. Notice again that a —>■ as r —>■ for > 0. The function is 
always zero at the initial singularity (a = 0) and increases towards its maximum value <fi* 
where uj — ► 00. Because r remains bounded (r < 1) in an open universe, <fi will never reach 
(jj* if r* > 1. 

For t > r*, <p then decreases towards zero (which it attains for k > as r — > 00). This 
demonstrates that, although — > as — > 0* and c<j — > 00, this is not the late time attractor 
solution. Instead we require that the function Y must diverge as r — ► 00 and thus — > 0. 

Vacuum solution with Brans-Dicke and G.R. limits 
Consider the function u(<f>) such that 

2o;(0) + 3 = (2uj + 3) - y - . - . (3.40) 
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Clearly u — ^constant as — > oo, but is divergent at 
Considering only > 0*, initially, we have 



r W = ,^I±iln , (3.41) 



and thus 




* rr + r* 1 + A;r 2 



(3.42) 
(3.43) 



Because Y(<j>) is divergent at both = 0* and — > oo we again have two distinct branches 
according to whether 0/0 is increasing or decreasing. In the former case, when (3$ > 0, 
we find a slow branch where the initial general relativistic behaviour ~ 0* turns into the 
Brans-Dicke solution oc as r — > oo. For (3q < we have decreasing 0, the fast branch, 
and the behaviour is reversed as r — > oo. Note that the r — > oo limit is only achieved for 
k > and that the late time behaviour in open models always corresponds to the general 
relativistic behaviour with — > 0*(1 + t~p°) as r — > 1. 

For < 0* we see that 2c<j + 3 may reach zero. This allows far more complex behaviour. 

= 0*(l-(f) /3 °) , (3-44) 

a 2 = 41 1— . (3.45) 

0, r A) _ r /3o 1 + A;r 2 v 7 

Notice now that reaches zero when r = r*. This corresponds to the divergence of the scale 
factor a at a finite proper time. Thus, for instance, the closed model does not recollapse. 

C. Stiff fluid solutions 

We consider in this section the case where matter is described by a barotropic equation 
of state with 7 = 2 (stiff matter) which as we have seen describes the long wavelength modes 
of a massless scalar field. In terms of the same variables X and Y the field equations become 

(X') 2 + AkX 2 - (Y 1 Xf = AMcj) , (3.46) 



X" + AkX = . (3.48) 

The last equation is identical to the corresponding equation for the vacuum case, and thus 
X(i]) is given by the same expressions (Eq. (|3.21|) ). This also yields the first integral 
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X' 2 + AkX 2 = A , (3.49) 
which upon insertion into the first of the field equations leads to 



Y'X = ±\JA - 4M0 . (3.50) 

This requires that <fi < A /AM. Notice also that unlike the vacuum case A could be negative, 
but only if k = — 1 and is also negative. This gives one extra solution for X(j]) when 
k = — 1 in addition to the vacuum solutions where A = —A 2 , 

A 

X{rj) = --cosh27? , (3.51) 



or in terms of the variable r defined in Eq. (|3.22 ) 



1 - kr 2 

A ' W = ~ A T+kP- ' (3 ' 52) 
For k — +1 note that this corresponds to X oc cos 2t] which is equivalent simply to a 
different choice of the zero-point of r), but in the open model we have a qualitatively different 
behaviour when A < 0. X = a 2 <ft remains non-zero at all times and thus we can obtain 
non-singular models where a remains non-zero. 
Now, from Eq. ( |3.50|) we derive 



Y = I ^A - 4M0^ , (3.53) 



and thus it is useful to define 



Z Ws/ ^_^ = = ± /^. (3.54) 



where the right-hand side of this equation is just ±lnr as for the vacuum case. Thus, just 
as in the vacuum case we required the function Y(4>) to diverge as r — > or — >• oo, in the 
stiff fluid case we require Z(<f>) to diverge in these limits. We see that if 

2a; vac (0) + 3 _ 2cu(0) + 3 
A 2 ~ A- 4M0 ' 

the vacuum solutions for a(t) and <fi(t) with uj vslc (4>) carry over to the stiff fluid solutions 
for uj((p). The reason for the equivalence becomes more apparent when we discuss the 
conformally transformed picture in the next section. When A < we see that for 2a; + 3 > 
we find the vacuum equivalent 2u vac + 3 < which is why we did not find the non-singular 
open models in the vacuum case. 

The condition for a = now becomes 

6(kX 2 - M<(>) 

U = A - AkX 2 ' (3 ' 56) 

confirming that u > is compatible with a turning point for k < when A < 0. For k = 
where we must have A > 0, or as X — > 0, the condition becomes uj = —6M<f)/A. Thus the 
sign of cf) becomes crucial. As one might expect, if the gravitational mass, M/<j) is negative, 
the initial singularity can be avoided even for u > 0, while for M/<ft > the presence of the 
stiff fluid requires an increasingly negative value of u to avoid the singularity. 
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Stiff fluid solution in Brans-Dicke gravity 



Proceeding as for the vacuum case, we start by considering the lu = ivq =constant case 
which enables us to compare our results with the k = solutions existing in the literature. 
For A = +A 2 > AM<p we have 



A - VA 2 - AMcj) 

A + y/A^tM^ 



± In r + constant . (3.57) 



Notice that this is exactly the same result as found in the vacuum case with 2c<j vac (0) + 3 = 
(2a;* + 3)0*/(0* - 0) if we write 2cu* + 3 = (2u + 3) and 0* = A 2 /(AM). This is a 
demonstration of the equivalence between different vacuum and stiff fluid solutions given in 
Eq. (ET55D . 



Thus, for A > 0, 



A 2 T Jjf_ 
~M 



2 , (3-58) 

(3.59) 



~ A riV l + kr 2 ' 

where r* is the constant of integration chosen to coincide with the value of r for which 
reaches its maximum possible value 0* = A 2 /AM. For r > r*, decreases back towards 
zero. [See Fig. ||. 

If k — — 1, r is bounded and will never attain 0* if r* > 1. In this case remains a 
monotonically increasing function of r approaching 40*/ (1 + rf ) 2 as r — >• 1 and thus t tends 
to infinity. 

If on the other hand we consider A < 0, we find a solution for <b < 0: 



A 2 T?T? 



n — n 



M (rf - 



2 ' 
2 



(3.60) 



It is possible to see that these solutions corresponds to the u > —3/2 solution derived 
by Gurevich et al |L6| (after the necessary translation to their time variable; Gurevich et al 
use £ such that d£ = drj/a 2 ). Notice that a = at both r = and r = r*, demonstrating 
that a turning point can indeed occur for u > even in open or flat models in the presence 
of the stiff fluid if < 0. As r approaches r* from below the solution approaches Nariai's 



power law solution |1J], but it is clear that this is not the late time behaviour suggested by 
Gurevich et al but rather a recollapse at a finite proper time. The correct late time behaviour 
for expanding k = models is where they approach the vacuum solution as r — > oo with 
positive or negative. 

When A = —A 2 , possible only for k < 0, we have X(rj) given by Eq. ( |3.52| ) and 



, s l2u + 3 i /V-4M0- A 2 \ . 
A x Z(0) = 2 W tan" 1 -f = ±2 tan" 1 r + constant . (3.62) 
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This gives 

A 2 



sec 2 



M/ (c + ^tan-V) . (3.63) 
Thus < — A 2 /AM as required. However — > — oo whenever r = tan((7r/2 — c)/f3) leading 



to a = \JX/(j) — > 0. This can always occur when the arbitrary constant c is sufficiently close 
to 7r/2 regardless of the sign of k. 

Stiff fluid solution with Brans-Dicke and G.R. limits 

If we consider again the function 2a>(0) + 3 = {2oj$ + 3)0 2 /(0 — 0*) 2 this time in the 
presence of a stiff fluid, we can integrate Eq. ( p.54|) for > 0* to give 



m = p^l ^=^= in 



Vv4 2 - 4M0 - V^ 2 - 4M0* 



V^ 2 - 4M0 + V^ 2 - 4M0* 



(3.64) 



Here must be constrained to lie within 0* < < A 2 /4M and so can never reach the 
asymptotic Brans-Dicke limit as — > oo. We find 

(rf-r B ) 2 0, + rfr B (A 2 /M) 



where we have written 



(3.65) 



IA'-Ut*._3_ 
V A 2 2a>o + 3 V ; 

Thus = 0* at r = 0, and reaches a maximum of = A 2 /4M when r = t* (possible only 
for r* < 1 in the open model). At late times for k > 0, as r — ► oo, returns to the general 
relativistic result, — > 0*, — > 0. 

Once again for < 0* we find a considerably more complicated behaviour where we may 
have 0^0 for non-zero X. 



D. Radiation solutions 

The other case in which the non-vacuum equations of motion simplify considerably is 
where the energy-momentum tensor is traceless (7 = 4/3), i.e. a radiation fluid correspond- 
ing to the short wavelength modes of a massless field. As this case has been discussed 
elsewhere we will describe the behaviour only briefly for comparison with the stiff fluid 
case, while presenting our results in a more compact form in terms of the time coordinate 
r(r]). 

The field equations in the presence of radiation with density p = 3r/87ra 4 where T is a 
constant, become 

(X') 2 + 4A;X 2 - (Y'X) 2 = ATX , (3.67) 
(Y'X)' = , (3.68) 
X" + 4A;X = 2r . (3.69) 
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The final equation can again be integrated directly to give the first equation where (Y'X) 2 = 
A 2 =constant. Notice that unlike the stiff fluid case this constant cannot be negative. The 
general solution of the equation of motion for X is then 



X 



t(A + Tt) 

1 + kT 2 



(3.70) 



in terms of the time coordinate r(r/) introduced in Eq. Q3.22 ). 

The Brans-Dicke field is not driven by matter and we have the same integral for Y( 
as in the vacuum case, although we have a different X(rj): 



Y(cj>) = ± 



Adr] 



±ln 



Tt 



A + Tt 



constant . 



(3.71) 



The evolution of 4>(r)) is thus the same as for the vacuum case if we replace the function 
t(v) by 



Tt( V ) 



A + Tt(tj) 



(3.72) 



Note that in spatially flat or closed models as r — ► oo we find s — > 1, i. e. (j) approaches a 
fixed value. In open models as r — > 1 we have s — > T/(A + T). On the other hand at early 
times these solutions approach the vacuum solutions as s ~ (T/A)t amounts simply to a 
rescaling of the conformal time or, equivalently, the scale factor. 

It is now straightforward to write down the radiation solutions for the particular choices 
of u>(4>) given in the vacuum and stiff fluid cases. 

The value of the u at any turning point is now given by 



6{kX 2 - TX) 
A 2 - A{kX 2 - TX) 
6TX 



A 2 + ATX J 



for k = or X -> 



(3.73) 
(3.74) 



The denominator must always be positive (by Eq. ( |3.68| )) and thus we find again that to 
obtain a turning point with u > requires either k > 0, which corresponds to the usual 
recollapse in closed models, or X (and thus <fi) negative. 



1. Radiation solution in Brans-Dicke gravity 



,±/3 



(3.75) 
2 1 s^t{A + Tt) 

a = J, 1 + kT 2 • (3 ' 76) 

As is well known, this Brans-Dicke solution approaches the general relativistic solution with 
constant u at late times during the radiation dominated era. 
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2. Radiation solution with Brans-Dicke and G.R. limits 



When 2u(<j>) + 3 = (2uo + 3)0 2 /(0 - 0*) 2 we find 




(3.77) 



1 sl« t(A + Tt) 



a 2 = — i — . (3.78) 

0, tg> + s fh 1 + kr 2 K J 

Once again approaches a constant at late times however, unlike the stiff fluid case con- 
sidered earlier, this constant value may not be close to 0* so this need not coincide with 
U) — > oo. 



IV. CONFORMALLY TRANSFORMED FRAME 

It has long been realised that a theory with varying gravitational coupling such as scalar- 
tensor gravity must be equivalent to one in which the gravitational coupling is constant but 
masses and lengths vary Mathematically this equivalence can be shown by using a 

conformally rescaled metric 

/ fk \ 

(4.1) 




0o is just an arbitrary constant introduced to keep the conformal factor dimensionless. 
Written in terms of this new metric and its scalar curvature, R, the scalar-tensor action 
given in Eq. (|2.1|) becomes 



S = — / d x\ -g 

l07T JM 



- 167T ! -Ify^tb + f^rl ^matter 



2 

where we introduce a new scalar field ip(4>) defined by 




(4.2) 



\ ^— - • (1.3) 

The gravitational Lagrangian is reduced simply to the Einstein-Hilbert Lagrangian of gen- 
eral relativity, albeit at the expense of altering the matter Lagrangian. Thus we shall refer 
to this as the Einstein frame. 

The arbitrary dimensional constant 0o plays the role of Newton's constant, G = 05" 1 . 
In order to avoid changing the signature, the conformal factor relating the metrics must be 
positive. So for < we must pick 0o < giving a negative gravitational constant in the 
Einstein frame. Not surprisingly then, the usual singularitiy theorems need not apply even 
in the Einstein frame for < 0. Similarly, in the definition of ip we require 0o(2u; + 3) > 0. 
If this were not the case we could instead define a scalar field 



jT / , 2a; + 3 d<j> A . 

but this would have a negative kinetic energy density, again invalidating the usual singularity 
theorems by breaking the dominant energy condition. However for > and uj > —3/2 the 
FRW models must contain singularities in the conformal Einstein frame where a — » 0. 
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A. Vacuum solutions 



The field equations are then, at least in vacuum (£ ma tter = 0), just the usual Einstein 
field equations of general relativity plus a massless scalar field, if). In particular, in a FRW 
universe (which remains homogeneous and isotropic under the homogeneous transformation) 
we have 

~ 9 87T p k , . 

6 (po er 

d 2 ib ~ dib , 

£ + & = (4.7) 

where the scale factor in the conformal frame a = {4>/4> ) 1/2 a, dt = (0/0 o ) 1/2 ^ and H = 
(dd/dt)/d. (Note that t is the time in the conformal frame and not to be confused with the 
conformally invariant time, r], used earlier.) The massless scalar field behaves, as it must, 
as a stiff fluid with density p = p = (dip / dt) 2 /2. 

Notice now that the variables X and Y introduced in the previous section correspond to 
the square of the conformal scale factor and the scalar field ip respectively. 

(4.8) 




The equations of motion for the conformal scale factor written in terms of X and for if) 
written in terms of Y and derivatives with respect to the conformal time rj are then precisely 
Eqs. Q3T7| - |3T9D solved in Section PTTBj . 



We can solve explicitly for X and Y as functions of rj because the stiff fluid continuity 
equation can be integrated directly (as for any perfect barotropic fluid) to give p oc a~ 6 . 
These results are independent of the form of u((f>). A particular choice of u)(4>) determines 
how <p is related to the stiff fluid field ip. To obtain <f>(rf) we must be able to perform 
the integral in Eq. ( |4.9|) , and thus we also obtain the scale factor in the original frame, 
a = (V0) 1/2 5. 

B. Non vacuum solutions 

If we include the matter lagrangian for a perfect fluid in the original scalar-tensor frame, 
then there is a non-trivial interaction between this matter and the scalar field, if>, in the 
Einstein frame. 

-^matter = J -^matter • (4-10) 

Thus the matter energy-momentum tensor, defined in the Einstein metric, 
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-9 dg a b 

is no longer independently conserved, 



f ab = ( V-^matter ) , (4.11) 




V a T ab = , (4-12) 



unless it is traceless, i.e. vacuum or radiation. The conformally transformed density, p = 
(4>o/(p) 2 p, and pressure, p = (0o/0) 2 P, of the fluid become dependent on and thus -0, so 
while the fluid retains the same barotropic equation of state it is no longer a perfect fluid in 
general. Note however that the overall energy-momentum tensor of the matter plus the ip 
field must be conserved as guaranteed in general relativity by the Ricci identity. 
We have the usual general relativistic equations of motion in a FRW model 



(413) 

3 0o a 



M +I j 2 = _4 1 p + p + 3(p + p) 
di 3 0o 

and the interaction leads to a transfer of energy between the original fluid and the stiff (ip) 
fluid: 



S = -3H(p + p) + -i= J-^-(3p - p)S , (4.15) 



Again we find two cases in which the problem simplifies. Firstly for radiation (p = 3p) 
there is no interaction and both continuity equations can be directly integrated and the 
conformal picture contains two non-interacting fluids: 

(4.17) 
(4.18) 

This is precisely the case considered recently by Barrow |22| (although without explicitly 
invoking the conformal frame) and discussed in Section D [D|. 



8tt a 


A 2 


30o P 


456 


8tt „ 


r 


o 7 Prad 
300 


' 



The second case in which we can find exact solutions is where the original fluid is itself a 
stiff fluid (or massless scalar field) in which case although there is an interaction between the 
two fluids, their combined dynamical effect is that of a single perfect stiff fluidQ, or massless 
scalar field x sa Y : 



3 The combined energy-momentum tensor of two interacting fluids is equivalent to that of a 
single perfect fluid provided their velocity fields are parallel. This must be true if both fluids are 
homogeneous as is the case here. Futhermore as they are both stiff fluids, p = p, in this case, their 
total pressure must be equal to their total density. 
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This is why we find exactly the same equation of motion for the scale factor in the conformal 
frame, a 2 oc X, in the stiff fluid case as in the vacuum case. Notice now that in the conformal 
frame we must have A = +A 2 > for a positive energy density. The non-singular solutions 
found when k < and A < with a stiff fluid in the Jordan frame correspond to solutions 
with negative enrgy density in the Einstein frame. 

The continuity equation for the original fluid can always be integrated to give 



o 




(4-3 7 )/2 

(4.20) 



and so in the stiff fluid case we have 



We have a 2 oc X as a function of 77 and we must now perform the integral in Eq. (|3.54j) to 
obtain (j)(rj). The change in the relation between and the total stiff fluid density in the 
Einstein frame compared with the vacuum case is equivalent to a different choice of u(<p) 
(which relates to ip) as demonstrated in Eq. ( |3.55| ). The vacuum case can of course be seen 
as a special case amongst the stiff fluid solutions, where M = 0, and thus uj(<fi) = c<j vac (0). 

We can also obtain exact solutions for radiation and stiff fluid in the original Jordan 
frame as the radiation remains decoupled in the Einstein frame and the interaction is solely 
between the two stiff fluids in that frame. Thus the equation of motion for the scale factor 
in the conformal frame is exactly the same as in the radiation only case, Eq. ( |3.69| ), while 
the equation for is the same as in the stiff fluid case, Eq. ( |3.54|) . 



Stiff fluid plus radiation in Brans-Dicke gravity 



To solve for the evolution of Brans-Dicke models (where u =constant) in the presence 
of both radiation and a stiff fluid, the conformal frame is particularly useful as the evolution 
of the conformal scale factor, or X = a 2 /(pQ, is exactly the same as for radiation only 
(Eq. ( |3.69| )). The evolution of then follows directly from Eq. ( [4.21 ) as 



P 



so that, from the definition of ip in Eq. 




(4.22) 



'2^0 + 3 



In 



A - VA 2 - 4M0 

A+ y/W^tM^ 



= ±ln 


Ft 


A + Yt 





+ constant . 



(4.23) 



Rewriting this to give and thus a = yX/<j> yields 
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M (sf + s^f t{A + Tt) 



A 2 



s* s 



1 + kr 2 



where 



s{r)) 



Tr 



A + Tt 



(4.24) 
(4.25) 

(4.26) 



is a constant of integration and (3 = y3/(2uj + 3). Thus the behaviour is very similar to 
that seen for stiff fluid only, except that the variable s takes the place of r. 

Unlike r, s — ^constant at late times for k = (where s — > 1) as well as for k < 
(where s — > T/(A + T)). Thus the Brans-Dicke field becomes frozen in at late times in the 
flat FRW universe, dominated by the friction due to Hubble expansion driven by radiation, 
just as it is in the open FRW model where the expansion becomes driven by the curvature. 
Only in the closed universe does the dynamical effect of the stiff fluid remain important. 
Note however that the radiation delays the recollapse which occurs at rj > tt/2. This means 
that r = tan(r/) becomes negative, but the solution is still well behaved as s > and both 
the conformal Einstein and Jordan (for u > 0) scale factors recollapse, X,a —>■ 0, when 
rj = n + tan _1 (— A/T), as s — > oo. 

Notice once again that the presence of a stiff fluid in the Jordan frame just leads to 
solutions which would be obtained in the absence of the stiff fluid but with the modified 
^vac(0) given in Eq. (|3.55|). 



V. CONCLUSIONS 



We have shown how to extend the procedure recently proposed by Barrow [22] to obtain 
the solutions for general cu(0) scalar-tensor gravity theories with a stiff fluid in addition to 
vacuum or radiation solutions in a FRW metric. These two non-vacuum cases correspond to 
the extreme long and short wavelength modes respectively of a minimally coupled massless 
scalar field. We show that these solutions can be obtained due to the particularly simple 
evolution of the corresponding scale factor in the conformally related Einstein frame which 
is independent of the form of u(4>). 

In the presence of a stiff fluid the scale factor evolves like a scalar-tensor cosmology with 
a modified u(<j)) — ► co> vac (0), as defined in Eq. ( |3.55| ). This is because introducing a new 
scalar field modifies the relation between (f) and the total energy density in the homogeneous 
scalar fluid. For example the Brans-Dicke model (u =constant) in the presence of a stiff 
fluid evolves like a vacuum model with uj va _ c ((p) oc 0*/(0* — (f>). This significantly modifies 
the evolution of the Brans-Dicke field leading to an upper bound of < 0*. 

We find that even for functions u(4>) that diverge at a finite value of 0, this need not 
be a stable late time attractor for k = models, in contrast to Damour and Nordvedt's 
rule [2(| that uj —>■ oo is a cosmological attractor. Instead (due to the absence of the 
damping effect of matter with p < p, required by Damour and Nordvedt's result) we find 
that the late (or early) time attractor in vacuum, as a — >• oo (or a — » 0), is associated with 
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the divergence of the function Y{4>) oc / \fluj + 3d0/0. In the presence of a stiff fluid the 
function Z(<f>) oc f(\J(2u + 3)/(A 2 — 4M0)<i0/0 must diverge as a — > or — > oo. 

The stiff fluid solutions are expected to be of primary importance as the scale factor 
a — > 0. When spatial curvature is negligible (k — 0), the condition necessary for a turning 
point, d = 0, in the stiff fluid cosmology is simply oo = —6M<f>/A where A and M are positive 
constants of integration. In vacuum this reduces to u = and the sign of is irrelevant. 
In the conformally related Einstein frame we have seen that the evolution is simply that for 
a stiff fluid irrespective of the form of uj(<f)) and thus singularities are always present here 
provided uj > —3/2. Only for uj < — 3/2 does the energy density of the stiff fluid in the 
Einstein frame become negative and so non-singular behaviour becomes possible. 
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FIGURES 

FIG. 1. The function X, defined in Eq. ( |3.20| ) plotted against conformal time, eta. The solid 
line represents k = 0, the dotted line k = +1 and the short dashed line k = —1 models. The long 



dashed line is the non-singular function for k = —1 given in Eq.(3.52). 



FIG. 2. Vacuum solutions for Brans-Dicke cosmologies showing Brans-Dicke field eft and scale 
factor a against proper time in the Jordan frame for the fast branch where <f) = at a = 0. Again 
the solid line represents k = 0, the dotted line k = +1 and the dashed line k = — 1 models. 

FIG. 3. Same as Fig. || but showing the slow branch. 

FIG. 4. Graph showing the exponents, q~(uj) (the fast branch) and q+(uj) (the slow branch), 
of power law expansion for k = vacuum Brans-Dicke cosmologies. 

FIG. 5. Vacuum solutions for scalar-tensor gravity theory with 2uj + 3 = 9(fi*/((ft* — (f>) showing 
Brans-Dicke field and scale factor against proper time in the Jordan frame. Note that to — > oo is 
not a late time attractor. This is identical to the evolution of a Brans-Dicke model in the presence 
of a stiff fluid. 
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